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We report on a search for the supersymmetric partner of the top quark (scalar top) decaying into 
a charm quark and a neutrahno in pp coUisions at y^ = 1.96 TeV. The data sample, collected by 
the CDF II detector at the Fermilab Tevatron, corresponds to an integrated luminosity of 2.6 fb~^. 
Candidate events are selected by requiring two or more jets and a large imbalance in the transverse 
momentum. To enhance the analysis sensitivity, at least one of the jets is required to be identified 
as originating from a charm quark using an algorithm specifically designed for this analysis. The 
selected events are in good agreement with standard model predictions. In the case of large mass 
splitting between the scalar top quark and the neutralino we exclude a scalar top quark mass below 
180 GeV/c^ at 95% confidence level. 

PACS numbers: 14.80.Ly 

I. INTRODUCTION 

The standard model (SM) of elementary particles and fundamental interactions, however successful, is still incom- 
plete since it does not explain the origin of electroweak symmetry breaking and does not give an answer to the gauge 
hierarchy problem jl|. A possible extension of the SM, supersymmetry (SUSY) [2], solves these problems by intro- 
ducing a symmetry that relates particles of different spin. R-parity [4I conserving SUSY models also provide a prime 
candidate for the dark matter in the universe 3], namely the stable lightest supersymmetric particle (LSP). In these 
models, the left-handed and right-handed quarks have scalar partners, respectively denoted as ql and qa, which can 



mix to form scalar quarks with mass eigenstates qi,2- Several models [4] predict that this mixing can be substantial 
for the scalar top (stop), yielding a stop mass eigenstate (ii) significantly lighter than other scalar quarks. If the 
ti is sufficiently light, it can be pair-produced copiously in proton-antiproton collisions at center-of-mass energy of 
^/s = 1.96 TeV at the Tevatron. A stop in the mass range accessible at the Tevatron is expected to decay, among other 
possibilities, into a charm quark and the lightest neutralino (xi), which is often assumed to be the LSP. Assuming 
R-parity conservation, stops are produced in pairs. We consider a scenario where m^ < mi, + mj^+ (where x^ is a 
chargino), and mf_^ < rriw + mf, + rrij^o . Under these conditions the decay ti — >■ cxi is dominant. Results from searches 
using this final state at the Tevatron have been previously reported in 151, |2] . 

In this Letter, we report the search for ti — > cxi decays in pp collision data from 2.6 fb^^ of integrated luminosity 
collected by the upgraded Collider Detector at Fermilab (CDF II) at the Tevatron. The final state contains two c 
jets from the hadronization of the c quarks and features an imbalance in transverse momentum ("missing transverse 
energy" or ^t [7]) from the two LSPs escaping detection. 

II. ANALYSIS DESCRIPTION 

CDF II is a multipurpose detector, described in detail elsewhere [g|. The charged-particle tracking system consists 
of silicon microstrip detectors and a cylindrical open-cell drift chamber, both of which are immersed in a 1.4 T 
solenoidal magnetic field coaxial with the beam line. The silicon detectors provide coverage in the pseudorapidity Q 
range \r]\ < 2 and are used to identify events with long-lived particle decays. The drift chamber surrounds the silicon 
detectors, has maximum tracking efficiency up to \ri\ = 1, and is used for charged particle momentum measurements. 
Segmented sampling calorimeters, which surround the tracking system, are arranged in a projective tower geometry 
and are used to measure the energy of interacting particles in the region I77I < 3.6. Muon candidates are identified by 
drift chambers, which extend up to I77I ~ 1.5, and are located outside the calorimeter volume. Jets are reconstructed 
from the energy depositions in the calorimeter cells using an iterative cone jet-clustering algorithm [9], with a cone 
size of radius R — ■\/(A0)2 + (Arj)'^ = 0.7 7]. Energy corrections [10] are applied to account for effects such as non- 
linear calorimeter response, underlying event, and the position of the interaction point, that influence the measured 
transverse jet energy. 

Candidate events used for this search are selected by an online event selection system (trigger) that requires 
l^T > 35 GeV and two jets. Further selections are applied off-line to remove accelerator-produced and detector- 
related backgrounds as well as cosmic-ray events. After off-line event reconstruction, the events are required to have 
l^T > 50 GeV, and at least two jets with \ri\ < 2.4 and Et > 25 GeV. The highest-E'^ jet is required to have 
Et > 35 GeV and at least one of the selected jets is required to be in the region \r]\ < 0.9. 

The hadrons in jets coming from & or c quark fragmentation, heavy-flavor (HF) jets, have a measurable flight path, 
yielding secondary vertices relative to the pp interaction point (primary vertex). We require the events to have exactly 
one jet identifled as a HF jet by the secondary- vertex tagging algorithm |11| . since requiring two HF-identified jets 
enhances the sample with events containing two b quarks and reduces the signal acceptance. This criteria defines 
a preselection which is used as the basis of background studies and further optimization for the signal sample, as 
described below. 

Dominant SM backgrounds are pair and single top-quark production, electroweak single and di-boson production, 
HF multijet production, and light-flavor multijet events where one of the jets is falsely tagged as a HF jet (mistag). 
The latter two background contributions are estimated from data. The ALPGEN |i2| event generator interfaced with 
the parton-shower model from the PYTHIA [1j| event generator is used to estimate the electroweak boson production, 
the MADGRAPH [IJ] generator is used to model the single top events, while the PYTHIA event generator is used to 
model the top-quark pair and diboson backgrounds. For the event generation the CTEQ5L (l5| parton distribution 
functions (PDFs) are used. Simulated events are passed through the GEANxS-based [la] CDF II detector simula- 
tion [17,] and are weighted by the probability that they would pass the t rigg er selection. The single top-quark and 
diboson contributions are normalized to the theoretical cross sections [18l420l ]. The contributions for the electroweak 
boson samples are normalized to the next-to-leading order cross sections calculated with mcfm [21| . We use the mea- 
sured top-quark pair production cross section of ati = 7.02 ± 0.63 pb [24]. Mistags are estimated based on the mistag 
rate |11| observed in inclusive jet data. The mistag rate is parametrized as a function of jet Et, ]??], secondary- vertex 
track-multiplicity, the number of primary vertices in the event, primary vertex z-position, and the scalar sum of Et 
of all jets in the event. To estimate the HF multijet background from data, we use a multijet tag rate estimator 
(MUTARE) described elsewhere [23i[2^]. The HF mufiijet prediction from MUTARE is scaled by a mufiiphcative 
factor that is obtained in a signal-free region. 

To avoid potential biases when searching for new physics, we test the various background contributions in distinct 
control regions that are defined a priori. The three control regions used to validate the SM prediction are denoted as 
multijet, lepton, and pre-optimization regions. The multijet control region is defined to have the second leading Et jet 



TABLE I: Comparison of the total number of expected and observed events in the control regions. The total uncertainty is 
computed by taking into account the (anti)correlations between the partial uncertainties. In the lepton region the HF multijet 
prediction is scaled to match the observed events in order to perform shape comparisons only. 



Regions: 


Multijet 


Lepton 


Pre-optimization 


WjZ + jets 


457 ± 190 


375 ± 156 


1551 ± 644 


Diboson 


17 ±2 


45 ±5 


118 ±13 


Top pair 


188 ± 21 


547 ± 60 


870 ± 96 


Single top 


11±2 


71 ±10 


130 ± 19 


HF multijets 


75407 ± 23376 


268 ± 83 


12935 ± 4010 


Light-flavor jets 


65839 ± 8427 


720 ± 92 


7741 ± 991 


Total expected 


141919 ± 24849 


2026 ± 208 


23345 ± 4182 


Observed 


143441 


2026 


22792 




50 100 150 200 250 300 350 400 

f T [GeV] 

FIG. 1: Distribution of Et in the pre-optimization region. SM prediction (stacked histograms) and observed distribution 
(dots) are shown, where HF multijets and light-flavor jets are predicted from data. 



direction (J2) aligned with the IfiT Q; where aligned means A(I)(^T:J2) < 0.4 rad. This HF multijet enriched region is 
used to obtain the MUTARE parameterization to predict the HF multijet background in the other control and signal 
regions. The lepton control region is defined to have J2 direction not aligned with the '^t i^4'(^T,J2) > 0.7 rad) 
and at least one isolated charged lepton (e or fi) with pr > 10 GeV/c . This lepton region is used to validate the 
modeling of the top quark and electroweak W/Z boson backgrounds. The pre-optimization control region is defined 
to have the leading and second leading Et jets not aligned with the '^t and to have no identified leptons, where the 
lepton definition is the same as in the lepton control region. We studied kinematic distributions such as jet Et, I^t, 
as well as the total numbers of events predicted. All these quantities are in agreement with observations in the three 
control regions. As an example, the }^t distributions for the pre-optimization region are shown in Fig. [T] Table U 
presents the expected event yields for various background processes in the control regions under study. 

We optimize the sensitivity to stop production by applying additional event selection criteria using a set of cuts 
suppressing the HF multijet background. We select events with exactly two jets, as expected from the signal, and 

fulfilling the condition A(f>(^ t ,^ t ) < ^/2- This variable is the angular difference between the calorimetry-based ^t 

and the same quantity calculated using tracks (^t ) |3|- When the I^t in the event is real, these two quantities 
are usually aligned in 0. However, when the I^Jt comes from calorimetry mis-measurements, as in HF multijet events 
with no real I^t, the angular difference between the two quantities is randomly distributed. 

To further improve the sensitivity we apply a neural network (NN) trained with the tmva package [25|, to reduce 
the remaining HF multijet background. We train the NN using jet Et, jet 77, the minimum A0 between the ^t and 

'T,lpT ) 



any of the selected jets, A(/)(7 



A(j){ji,J2), l^T, I^T^' ^^'i th^ summed Et of all the jets in the event. We 



TABLE II: Number of expected and observed events in the signal region before and after CHAOS application. Prediction for 
the signal point (m(ti) — 125 GeV/c^ , m(xi) = 70 GeY/c^ ) is also shown. Correlated and uncorrelated uncertainties in the 
total background and expected signal were treated separately in the analysis although they are combined here. 



Signal region Before CHAOS Final (after CHAOS) 

W/Z + jets 423.6 ± 185.0 60.9 ± 26.6 

Diboson 36.9 ± 6.5 10.7 ±1.9 

Top pair 61.9 ± 15.5 4.6 ± 1.3 

Single top 39.0 ± 9.7 3.2 ± 0.8 

HF Multijets 279.6 ± 208.3 20.4 ± 15.2 

Light-flavor jets 658.3 ± 259.6 32.2 ± 12.7 



Total expected 


1499.3 ±277.1 


132.0 ±24.4 


Observed 


1496 


115 


ti signal 


250.0 ±66.2 


90.2 ±23.9 



choose a reference signal point with m{ti) — 125 GeV/c^ and to(xi) — 70 GeV/c^ to perform the optimization. The 
signal acceptance is obtained using the pythia event generator and CTEQ5L PDFs. Total signal yields are normalized 
to the NLO production cross section determined with the prospino event generator [2^ and the CTEQ6M [13, HI] 
PDFs. The uncertainty of the NLO production cross section is estimated to be 20%, arising from the scale dependence 
and the uncertainties on the PDFs. The NN output lies within —1 and 1, where the background peaks at —1 and the 
signal peaks at 1. We define our signal region as events with NN output scores > 0. 

The final stage in the optimization is the application of a charm hadron analysis oriented separator (CHAOS) 
technique [2J], explicitly designed for this analysis to obtain a sample enriched in c jets. CHAOS is a NN producing 
a two-dimensional output and trained with the SNNS v4.3 package [23] to determine whether a jet identified as HF 
has been produced from the hadronization process of a light quark falsely tagged as a HF jet, a b quark, or a c quark. 
The two-dimensional output structure allows the separation of the three different targets during the same training 
process. Depending on the flavor of the original parton, the jet identified as HF and its secondary vertex have different 
characteristics, mainly related to the tracking. Using properties of the tracks forming the secondary vertex and the 
tracks of the jets within a neural network, CHAOS allows enhancement of the jet selection with a desired fiavor, 
in particular c jets. We apply CHAOS to the jet identified as HF and we find that the optimal cut has a selection 
efficiency [24] of 34% for c jets, 7.3% for b jets, and 4.9% for light jets. 

The systematic uncertainties on the signal and the background predictions, taking into account correlated and 
uncorrelated uncertainties, are studied. Correlated uncertainties, affecting both the background prediction and signal 
acceptance, are dominated by the uncertainties on the performance of the 6-tagging algorithm and CHAOS, which 
are 4.4% [llj and 9.2% respectively, and the luminosity (6%) Q. Uncorrelated systematic uncertainties on the 
background predictions are dominated by uncertainties on the MUTARE parameterization (30%), the mistag rate 
(16% [Tl!] for light-fiavor multijets), the top-quark pair-production cross section (11%), the single top-quark production 
cross section (13%), and the diboson production cross section (10% for WW/WZ and 20% for ZZ). The uncertainty 
on the normalization of the boson plus HF jets to the total inclusive boson production cross section translates into 
a 10% uncertainty in the SM predictions. Correlated and uncorrelated uncertainties are evaluated separately and 
combined in quadrature. 

III. RESULTS 

The signal region is analyzed after the background predictions are determined. We observe 115 events, where 
132.0 ± 24.4 are expected from background, as summarized in Table IhI where yields before selecting events based on 
CHAOS output are also shown. Since no significant deviation from the SM prediction is observed, the results are 
used to calculate a 95% C.L. exclusion limit for the ti pair production cross section. 

We have used the differences in shape of the NN output to set the limits. These limits are computed using a 
Bayesian likelihood method fSlj with a fiat prior probability for the signal cross section and Gaussian priors for the 
uncertainties on acceptance and backgrounds. Figure [2] shows the expected and observed limits as a function of m(ti) 
for a neutralino mass of 80 GeV/c^ . 

We exclude, assuming BR(ii — )• ex?) = 100%, ti masses up to 180 GeV/c^ at 95% C.L. In addition, a 95% C.L. 
limit is obtained in the mass parameter plane of the model. Figure [3] shows the excluded region in the stop-ncutralino 
mass plane of the analysis, compared with results from previous analyses |5l. la. l30l]. The limit obtained with the 
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FIG. 2: Observed (solid line) and expected (solid line with dashed band) 95% C.L. upper limit on the stop cross section (solid 
line with shaded band) as a function of the stop mass for an assumed value of the neutralino mass. The shaded band denotes 
the uncertainty on the NLO stop pair-production cross section. 
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FIG. 3: Excluded region at 95% C.L. in the m{ii)-m{xi) plane. The result is compared to the previous results from CDF [3], 
from DO [^, and from LEP ^ experiments at CERN. 

present analysis improves the results of previous searches using a similar topology, and represents the world's best 
limit in the region of large mass splitting. 



IV. SUMMARY 



To summarize, we have searched for the pair production of stop decaying into a charm quark and a neutralino, 
in 2.6 fb~^ of CDF Run II data. We observe 115 candidate events, which are in agreement with SM background 
expectations of 132.0 ± 24.4 events. No evidence for stop is observed, and we exclude a region in the stop and 
neutralino mass plane at 95% C.L. as shown in Fig. [3] Assuming BR{ii — )• ex?) = 100%, we exclude stop masses up 
to 180 GeV/c^ at 95% C.L. for a neutralino mass of 90 GeV/c^ . 



Acknowledgements 

We thank the Fcrmilab staff and the technical staffs of the participating institutions for their vital contributions. 
This work was supported by the U.S. Department of Energy and National Science Foundation; the Italian Istituto 
Nazionale di Fisica Nucleare; the Ministry of Education, Culture, Sports, Science and Technology of Japan; the 
Natural Sciences and Engineering Research Council of Canada; the National Science Council of the Republic of 
China; the Swiss National Science Foundation; the A. P. Sloan Foundation; the Bundesministerium fiir Bildung und 
Forschung, Germany; the Korean World Class University Program, the National Research Foundation of Korea; the 
Science and Technology Facilities Council and the Royal Society, UK; the Russian Foundation for Basic Research; 
the Ministerio de Ciencia e Innovacion, and Programa Consolider-Ingenio 2010, Spain; the Slovak R&D Agency; the 
Academy of Finland; and the Australian Research Council (ARC). 



[8] 
[9] 

[lo; 
[11 

[12 

[13 

[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22' 
[23^ 
[24' 
[25^ 
[26; 

[27 
[28' 
[29 
[30 
[31 



[1] E. Gildener, Phys. Rev. D 14, 1667 (19 76) 

[2] For a review of SUSY, see S. P. Martin, |hep-ph/9709356 and references therein 



[3] H. Goldberg, Phys. Rev. Lett. 50, 1419 (1983); J. Ellis et ai, Nucl. Phys. B238 (1984) 453. 

[4] A. Bartl, W. Majerotto, and W. Porod, Z. Phys. C 64, 499 (1994) [Erratum-ibid. C 68, 518 (1995)]. 

[5] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 76, 072010 (2007). 

[6] V.M. Abazov et al. (DO CoUaboration), Phys. Lett. B 665, 1 (2008). 

[7] We use a cylindrical coordinate system with its origin at the center of the detector, where d and (j> are the polar and azimuthal 
angles, respectively, and pseudorapidity is r; = — ln(tan(|)). The missing Et (Et) is defined by 'Et ~ — '^^ Et^i, where 
i = calorimeter tower number and lii is a unit vector perpendicular to the beam axis and pointing at the i**^ calorimeter 
tower. 'Et is corrected for high-energy muons and jet energy. We define Et = [Et]- The I^t is defined as the negative 
vector sum of track pr's requiring tracks with pt > 0.5 GeV . We define E^^ — \Eq 



'T 



D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 032001 (2005). 

G. Arnison et al. (UAl CoUaboration), Phys. Lett. B 123, 115 (1983); A. Bhatti et al, Nucl. lustrum. Methods A 566, 

375 (2006). 

A. Bhatti et al. (CDF Collaboration), Nucl. lustrum. Methods A 566, 375 (2006). 

D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 052003 (2005). 

M. L. Mangano et al, J. High Energy Phys. 07 (2003) 001. We use ALPGEN v. 2.10. 

T. Sjostrand et al, Comput. Phys. Commun. 135, 238 (2001). We use pythia v. 6.216. 

T. Stelzer and W.F. Long, Comput. Phys. Commun. 81, 357 (1994). 

H. L. Lai et al, Eur. Phys. J. C 12, 375 (2000). 

R. Brun et al. Tech. Rep. CERN-DD/EE/84.1, 1987. 

E. Gerchtein and M. Paulini, eConf C0303241, TUMT005 (2003). 

B. W. Harris et al, Phys. Rev. D 66, 054024 (2002); Z. Sullivan, Phys. Rev. D 70, 114012 (2004). 

M. Cacciari et al, J. High Energy Phys. 04 (2004) 068. 

J. R. CampbeU and R. K. Ellis, Phys. Rev. D 60, 113006 (1999). 

J. CampbeU and R. K. Ellis, MCFM- Monte Carlo for FeMtobarn processes, |http://mcfm.fnal.gov/[ We use mcfm v. 5.8. 

A. Lister (CDF and DO Collaborations), arXiv:0810.3350 

T. Aahonen et al (CDF Collaboration), Phys. Rev. Lett. 102, 221801(2009). 

M. Vidal, Ph.D thesis, C IEMAT, FERMILAB-THESIS-2010-08, 2010. 

A. Hocker et a?., |arXiv:p hysics/0703039 



W. Beenakker et al, Nucl. Phys. B492, 51 (1997); W. Beenakker, R. Hopker, and M. Spira, ' hep-ph/9611232[ We use 

PROSPINO V. 1.0. 

J. Pumplin et al, J. High Energy Phys. 07 (2002) 012. 

D. Stump et al, J. High Energy Phys. 10 (2003) 046. 



A. Zell et al, http://www.ra.cs.unituebingen.de/SNNS/contact.html 
LEPSUSYWG/04-01 .1,. http://lepsusyweb.cern.ch/lepsusy/, 



J. Heinrich et al, arXiv:physics/0409129 



